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SMILANSKY’S MODEL OF IRREVERSIBLE QUANTUM 
GRAPHS, I: THE ABSOLUTELY CONTINUOUS SPECTRUM 

W.D. EVANS AND M. SOLOMYAK 


Abstract. In the model suggested by Smilansky [7| one studies an opera¬ 
tor describing the interaction between a quantum graph and a system of K 
one-dimensional oscillators attached at several different points in the graph. 
The present paper is the first one in which the case AT > 1 is investigated. 
For the sake of simplicity we consider AT = 2, but our argument is of a 
general character. In this first of two papers on the problem, we describe 
the absolutely continuous spectrum. Our approach is based upon scattering 
theory. 


int 


1. Introduction 

In the paper [7] U. Smilansky suggested a mathematical model to which he 
gave the name “Irreversible quantum graph”. In this model one studies the 
interaction between a quantum graph and a hnite system of one-dimensional 
oscillators attached at several different points in the graph. Recall that the 
term “quantum graph” usually stands for a metric graph T equipped with a 
self-adjoint differential operator acting on L^(r); see the survey paper jl] and 
references therein. In our case this operator will be the Laplacian —A. 

In Smilansky’s model one initially deals with two independent dynamical 
systems. One of the systems acts in L^(r) and its Hamiltonian is the Laplacian. 
Another system acts in the space K > 1 and is generated by the 

Hamiltonian Hose = Ylk=i^k where 


hk 



/c = 1,..., A; 


in [7j the oscillators are written in a slightly different form; one form reduces 
to another by scaling. In what follows the points in T are denoted by x and 
the points in by q = (gi,..., q^)- 
Consider now the operator 


1.1 


( 1 . 1 ) 


Ag — —A (K> / “b / ® Hose 
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in the space L^(r x It is defined by the differential expression 


1.2 


( 1 . 2 ) 


K 


AU = 


k=l 


d^U 

dqt"^ 


+ qlu) 


and is self-adjoint on the natural domain. The terms in (HID do not interact 
with each other. 

Interaction is introduced with the help of a system of “matching conditions” 
on the derivative f/j. at some points oi,... ,ox G T. One says that the fc-th 
oscillator is attached to the graph at the point Ok- The condition at the point 
Ok is 


1.4 


(1.3) [U'^]{ok,q) = akqkU{ok,q), k = l,...,K, 

where [/^](.) stands for the expression appearing in the Kirchhoff condition, 
well known in the theory of electric networks. When T = M (which is the only 
case we deal with in the main body of the paper), [/j,](.) is the jump of the 
derivative. 


1.5 


(1-4) [/'](o) = /'(o+)-/'(o-). 

The real parameter in dH expresses the strength of interaction between the 
quantum graph and the oscillator hk- The case «! = ... = ax = 0 corresponds 
to the operator Aq as in (HID. 

Sometimes we shall denote hy a, v the multi-dimensional parameters a = 
{«!,... v = {z/i,..., Let stand for the 

operator defined by the differential expression (j1.2^1 and the conditions ina. 
Usually, the values of Vk are fixed and we exclude them from the notation. On 
the other hand, we use the notation Ar.Q,;j, for this operator when it is necessary 
to reflect its dependence on the graph. 

The problem to be considered is the description of the spectrum of the 
dynamical system generated by the Hamiltonian Aa-^y. More specifically, it is 
to construct the self-adjoint realization of Aa-^y as an operator in the Hilbert 
space L^(r x M-^) and to describe its spectrum. 


Up until now, the problem has only been investigated for the simplest case 
K = 1. The first results were obtained in the paper [7j by Smilansky. Then a 
detailed study of the problem was carried out in the papers i, i and j3]. In 
PI. along with some new results, a detailed survey of the current state of the 
problem is given. 


On first sight, the problem might seem amenable to the perturbation the¬ 
ory of quadratic forms. Indeed, the spectrum ct(Ao) can be easily described 
by separation of variables and the perturbation in the quadratic form, which 
appears when passing from Aq to Aq, with a ^ 0, seems not to be too strong. 
However, this is not so: this perturbation turns out to be only form-bounded 
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but not form-compact, which makes it impossible to apply the standard tech¬ 
niques. So, the problem requires certain specihc tools which were developed in 
0 - iia and j3]. The most important of these tools is the systematic use of 
Jacobi matrices. 

It was found in the above mentioned papers on the one-oscillator problem 
that the character of the spectrum strongly depends on the size of a\ there 
exists some a* > 0 such that the absolutely continuous spectrum (Ja.cX-^a) 
coincides with (Ta,c.(Ao) if |a| < a* (in particular, it is absent if the graph is 
compact) and hlls the whole of M if |a| > a*. The dependence of the structure 
of the point spectrum ap{Aa) on a is also well understood. 

This is the hrst of two papers on the problem for K > 1 and in it we study 
the absolutely continuous spectrum; in our other paper pp the point spectrum 
is investigated. This division is natural, since the technical tools used in each 
part are different. We address the simplest situation, when T = R and K = 2, 
but our argument is of a rather general character and we hrmly believe that it 
applies to a wide class of graphs and to any K. However, in the general case, 
the calculations become more complicated and this obscures the main features 
of the argument. 

We hrst describe informally the main ideas lying behind our approach. 

The effect of adding one more oscillator to a system with K oscillators is 
twofold. Firstly, the total dimension of the set T x R^ increases by one which 
certainly affects the spectrum. Secondly, there is some effect coming from the 
additional matching condition (El at the point ok+i- This second effect disap¬ 
pears if we take ax+i = 0. Indeed, then the variable qx+i can be separated and 
the operator decomposes into the orthogonal sum of simpler operators. More 
exactly, denote by A the operator which corresponds to the conhguration with 
the {K + l)-th oscillator removed, 

A — A 

Then it is easy to see that 

__m ^ 


nGNo 

This orthogonal decomposition yields the complete description of the spectrum 
of the operator on the left-hand side, provided that the spectrum of A is known. 

The key observation which allows one to solve the general problem is that 
the interaction between the oscillators attached at different points is weak. For 
K = 2 this observation leads to the conclusion that the study of a{Aa^^a 2 ]ui,u 2 ) 
can be reduced to the same problem for the operators Aaifi-ui,u 2 Ao,a 2 ;,yi,i/ 2 . 
Due to the equality EH) this reduces the problem to the study of the spectra 


1.6 (1.5) 


A 




= } (A + 4.(u + 1/2)). 
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St. 1 


st2 


of two operators, and each corresponding to the case of only- 

one oscillator. Since the latter case is already well nnderstood, we obtain the 
desired resnlts for onr more complicated case. 

An accnrate realization of this idea is different for the point spectrnm and 
for the absolntely continnons spectrnm. In the present paper we concentrate 
on the absolntely continnons spectrnm. Here an important correction to the 
above scheme is necessary: the stndy of (ya.c.{,Aa^^aT,y\,v-^ does not rednce to 
the stndy of oa.c.^Aa^.v-^ and oa.c.{AaY,v^ for the same graph F. Rather, we 
have to divide F into two parts, F = FiUF 2 in snch a way that Oj G F^ and 
Oj ^ Fg.j. Then cra.c.(Ar;ai,a2;i^i,!^2) can be expressed in terms of (Ja.c.{AYy^ocy,v^, 
j = 1,2. The paper [T] is devoted to the stndy of the point spectrnm. There 
snch a partition of F is nnnecessary. 

We nse the following notation. We write No for the set {0,1,...}. The diag¬ 
onal operator in an appropriate f'^-space, with the diagonal elements Oq, Oi,..., 
is denoted by diag{a„}. We apply similar notation for the block-diagonal oper¬ 
ators. The notation {bn}) stands for the Jacobi matrix whose non-zero 

entries are jn,n = cin and jn,n+i = jn+i,n = bn- If A is a self-adjoint operator in 
a Hilbert space, then ct(A), cTa.c(A), crp(A) stand for its spectrnm, absolntely 
continnons (a.c.) spectrnm and point spectrnm respectively. We nse the sym¬ 
bol ma.c(/; A) for the mnltiplicity fnnction of the a.c. spectrnm. The symbol 
Gi stands for the trace class of compact operators. 

Other necessary notations are introdnced in the conrse of the presentation. 

2. Statement of the problem. Results 

2.1. The operator Aq. As was mentioned in the introdnction, we present 
onr argnment for the graph F = R and K = 2. We choose the points oi = 
1, 02 = —1 and denote the coordinates in R^ by g_ and the parameters by 
a = {q;+, O-}, V = {z/+, p_}. The Laplacian on F is jnst the operator —(P/dx^ 
with the Sobolev space iF^(R) as the operator domain. The operator A^^y acts 
in the Hilbert space = L^(R^) and is dehned by the differential expression 

2 2 

(2.1) AU = A,U = -U'i. + + qlU) + + qlU) 

and the matching conditions (cf. (I1.4|l l 

(2.2) [ U {]{± l , q +, q -) = a±g±t/(±l, g+, g_). 

So, in the notation of the introduction, we are dealing with the operator 
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sa. 0 


sa. 1 


sa. 2 


sa. 3 


opO 


addl 


However, as a rule we use the shortened notation Aq. Note that the replace¬ 
ment a± ^ —a± corresponds to the change of variables q± ^ —q± which does 
not affect the spectrum. For this reason, we discuss only a± > 0. 

The structure of the differential expression A makes it natural to decompose 
the function 17 in a double series in terms of the normalized Hermite functions 
An, namely 

(2.3) V{x,q+,q-)= ^ Mm,n(a:)Am(g+)An(g-), 

m,n€No 


which is hereafter represented by 17 ~ {um,n}- The mapping U h->• {um.n} is 
an isometry of the space L^(R^) onto the Hilbert space S) = L^(M)). We 

evidently have 7117 {L^,nUm,n} where 

(2.4) {Lra,nU){x) = -u”{x) A rm,nU{x), X 7 ^ ±1; 


(2.5) rm,n = vX{rn +1/2)+ v‘^_{n +1/2), m,n ^Nq. 

The conditions at x = ±1 reduce to 


( 2 . 6 ) 


«m,n](l) = ^ (Vm + lu^+i^nil) + VmM„^_l,n(l)) ; 


[^m,n] ( T ^Um,n+l{ 1) T 1)^ • 


To derive the conditions (ESI) from (O, one uses the recurrency equation for 
the functions 

Vn+ lAn+l(g) - V^qXniq) + y/nXn-liq) = 0. 


2 .2. Operator Ag. The operator Ag := Ao g;i,^_,^_ admits separation of vari¬ 
ables and we get 

© 

(2.7) Ag = y~~^(Hg -|- r^.n); 

m,n 

where Hg is the self-adjoint operator —d?/dx^ in L^(M). This leads to the 
complete description of the spectrum ct(A g), namely, that it is purely a.c. 
and hlls the half-line [rg g, cx)) = [(i/^ -|- v//)/2,oo). The expression for the 
multiplicity function ma.cX^j 2^o) immediately follows from ()2.7j) . but is omitted. 
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opa 


0 




sa.thm 


circ 


2.3. Domain of A^. It is convenient to describe the domain of the self-adjoint 
realization of the operator Aq, in terms of the decomposition (Q. Dehne the 
set Dq, as follows. 

Definition 2.1. An element U ~ {um,n} lies in !)„ if and only if 
1 - Um,n £ H^{R) for all m,n. 

2. For all m, n the restriction of Um,n to each interval (— cxo, —1), (—1,1), (1, cx)) 
lies in and moreover, 




I Lm,n^m,n \ 


^dx < oo. 


3. The conditions (EH) are satisfied. 

Along with the set Dq, define its snbset 


D* = {D e Dq : D ~ {um^n} finite} 

where by finite we mean that the seqnence has only a finite nnmber of non-zero 
components. Denote by A* the operator in defined by the system 

(Q on the domain D*, 

Lemma 2.2. The operator A* is symmetric in S). Its adjoint coincides with 
the operator Aq considered on the domain DqI 


(A^)* = Aq. 

The proof is a straightforward modification of that for (5.2) in j3]. 

Theorem 2.3. For any «+, a_ > 0 the operator Aq is self-adjoint. 

The proof is given in section ^ Theorem 12.31 and Lemma 12.21 show that Aq 
is the nniqne natural self-adjoint realization of the operator, defined by the 
differential expression m and the matching conditions ()2.2j) . 


2.4. Absolutely continuous spectrum of the operator Aq. Below we 
construct an operator A° whose a.c. spectrum admits a complete description. 
Then we show that the a.c. spectra of both operators Aq and A° coincide, 
including the multiplicities. 

As a first step, let us consider two operators, A+^ and A~_. The operator 
Aq_^, say, acts in the space x M^) and is defined by the differential 

expression (ED), the matching condition (O at the point oi = 1 and the 
Dirichlet condition U{0, g+, g_) = 0. The deffnition of A~_ is similar, with M_|_ 
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replaced by M_ and the point oi = 1 by 02 = — 1. By separation of variables, the 
operators can be identihed with the orthogonal sum of simpler operators: 


one. 2 


( 2 . 8 ) 


^a+ — ^ + Z/^(n + 1/2)), 

nGNo 

0 

A„_ = (Ar_.„_.,^_ + z/^(m + 1/2)). 

mgNo 


one. 3 


Hence, both operators are self-adjoint. The direct sum 

(2.9) Al = := Ay e A;_ 

is a self-adjoint operator in the original Hilbert space .^. 


The following theorem is the main result of the paper. Its formulation in¬ 
volves the notion of wave operator, which is one of the basic notions in math¬ 
ematical scattering theory; see e.g. |S|, M or m 


one. tl 


Theorem 2.4. For each of the pairs (Aq,, A°), (A°, Aq), there exist complete 
isometric wave operators. In particular, the absolutely continuous parts of Aq 
and A° are unitarily equivalent. 


Theorem 1231 and the formulae (O, (EU reduce the study of <Ja.c.{h^a) to 
the similar problem for the case of only one oscillator. The latter problem 
was solved in j3] and uni. The next statement collects, for the particular case 
we need, the results of section 3 in PI; see also Theorem 5.1 and remarks 
in section 9 of j3]. In both papers it was assumed that u = 1, and we arrive 
at the formulation below via scaling. By default, we take ma.c.(^;A) = 0 if 
i i o-a.c.(A). 


prop 


Proposition 2.5. (The case of one oscillator.) Let T = M_|_ and o 
T = R_ and o = —1. Then 
1) <Ta.c. (Ao;i.) = [pV2, 00 ); 

vcia.c.il] Ao;i/) = n for — j2 < I — i/‘^n < p^/ 2, n G N; 


2) ifO<a< u\/2, then 

I^a.c. (Aq,.^) <Ta.c.( Aq.j/) [r' /2,CX)), 


1, or 


Rla.c. (^j ILIqi.c. Ao;i/), 

3) if a = u\/2, then 

<Ta.c. ( Aq.;/) [ 0 , CX)), TtTa.c. (1) Aq,.^) tTla.c. (Ij Aqjj^) 1, VI ^ 0, 

4) if u\/2 < a < oo, then 

I^a.c. ( Aq;j/) M, Aq-j/) Rla.c. (^j Aq.j/) 1, V/ G M. 
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adds 


one.t2a 


kato 


Now we are in a position to present the final formula for the function 
iTia.cX^j thus for our Original operator A^. 


( 2 . 10 ) 


^a.c. ^ ^ ^a.c.(^ h'_(^n 1/2), A]u^.Q,^.jy^) 

nGNo 

+ + 1/2); 


mSNo 

Combining the equality (|^:tud with Proposition 12.51 we obtain the following 
description of the a.c. spectrum of the operator Aq, for any «+,«_ > 0. 


Theorem 2.6. Let A„ = Aa-,^ be the self-adjoint operator defined by the dif¬ 
ferential expression (ED) on the operator domain Dq. 

1) then 

cra.c.(A„) = [ro,o, oo) = [(z/+ + z/^)/2, cx)). 

2) Let Q!+/z^+ = and a^jv^ < \/2, or a^jv^ = a/2 and a^/u^ < \pl. 
Then 

aa.cfAa) = [iyl/2,oo) or aa.c. (A„) = [z///2, cx)) 

respectively. 

3) Let a+/i^+ = a_lv_ = a/2, then 

CXa.c.(A„) = [0, CX)). 

In all the cases 1-3 the multiplicity function tria.c//; A^,), given by the 
equality fITTUD . is finite for all I G cTq.c/Aq,). 

4) Let max{a+/i>+,a-/v-) > \pl. Then 

^a.c. (Aq) M, IT^a.c. (^) Aq) = CX). 

In connection with this theorem, we would like to emphasize that the ex¬ 
istence of the wave operators established in Theorem 12.41 gives much more 
information about the operator Aq than just the description of its a.c. spec¬ 
trum. 

For the proof of Theorem 12.41 we use the following classical result due to 
Kato, see Theorem 6.5.1 and Remark 6.5.2 in m 

Proposition 2.7. Let A, A° be self-adjoint operators in a Hilbert space. Sup¬ 
pose that for some natural number p the inclusion 

{A° - L)-P - {A° - L)-P e 6i 

is satisfied for all non-real i G C. Then the complete isometric wave operators 
exist for both pairs A, A° and A°, A. 

In our case the conditions of Proposition 12.71 turn out to be fulfilled with 
p = 3. This is the result of the following statement whose proof is our main 
technical goal in this paper. 
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one. t2 




goal 


Theorem 2.8. For any non-real i G C one has 

(2.11) (A°-£)-3-(A,-i)-3e6i. 


Theorem 12.41 is a direct consequence of Theorem 12.81 

The proof of Theorem 12.81 is rather long and requires some preparatory work. 


_, 3. Auxiliary material 

au I 

In this section we present some elementary technical material concerning the 
equations 


res.Ig 
res.Igo 


(3.1) 

(3.2) 


-u" + Cn = /. 
-k" + CT = 0. 


where C = 7 + is a complex parameter. We need this material for the proofs 
of both our main technical results, Theorem 12.31 and Theorem 12.81 

We assume that 7 > 0, and are mainly interested in estimates which are 
uniform with respect to C,. 

3.1. Homogeneous equation. Let be the two-dimensional space of func¬ 
tions on M which are continuous, vanish as |a:| —cxo, and for a: 7 ^ ±1 satisfy 
the equation (Q. We choose the following basis in 




0 , 

sinh <^(x+l) 
sinh 2C ’ 
g-C(x-l) 


X < —1, 

—1 < X < 1; 
X > 1; 


Then 


(x) = ‘fifi-x). 


[I] (3.3) ^yi) = ^y-i) = i, ^+(-i) = »’c (!) = ')• 

Just for this reason this basis is more convenient than the “natural” basis 
consisting of the functions Note also that 


[(V^C )'] (± 1 ) 


2C 

1 -e-A’ 




2 Ce-^^ 

1 — e“A 


and hence, for all n G 


E (3.4) 


M(l) 

M(-i) 


^-^(n(l)-e 

H-1) - e~'‘^v(l)} . 
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lac 


2c 

nhom 


res. 1 


res.lb 


A standard calculation shows that for the norm and scalar product in L" 


(3.5) llv^jr = ||(^cH =Tr + 


7 ^ sinh 47 — 5 ^ sin 45 


27 4(sinh^ 27 + sin^ 25) 
=0{e-‘^'^), 7^cx). 


= 7 ^ + o(e ^^), 


This shows that for 7 large the chosen basis is “almost orthogonal”. It follows 
that the two-sided estimate 


(3.6) Co^ 7 ||n|h < IC+h-b < Co 7 ||n|n v = C+ipJ + e 


with some Cq > 1 is satished uniformly in any half-plane 7 = Re^ > 7 o > 0 . 

Now we turn to the subspace formed by the functions v G satisfying 
an additional condition n(0) = 0. The functions 


{ 0, x < 0, 

iSf. o<x<i, 7'-w = 7'+(-i.) 

g-C(a:-l)^ X > 1; 

form a natural basis in For the functions the equalities (13. 3 j) are 
satished, and instead of we have 

(3.7) [fy](±l) = - ^_^^_^^ n(±l), VugT^T 

Similarly to (ESD, we hnd that 

= 7 "^+ 0 (e“ 2 ^), = 0 (e“'^), 7 ^ cx). 

As a consequence, we conclude that an analogue of (ESI), with the functions 
replaced by is valid for n G 

A straightforward calculation shows also that 
(3.8) -(pjll = 0(e“^), 7 ^cx). 


3.2. Non-homogeneous equation. Here we discuss the equation m with¬ 
out the matching conditions at x = ±1 or, equivalently, under the conditions 
of the type dZSD with 0 = 0. Then the solution is given by 

(3.9) = 

7r 


The solution of the same equation dUD subject to the condition n( 0 ) 


(3.10) 


M^(x) 


(2C) ^ /r+ (e - e x > 0; 

( 2 C)“^ f{t)dt, x < 0 . 


0 is 
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res.If 


add 


The difference — U(^ is given by a rank one operator, 

(3.11) ul{x)-u^{x) =[ f{t)g(;{t)dt, g^{x) = 

Jr 

Note that = 7 “^. 

3.3. Dependence on the additional parameters. We are particularly in¬ 
terested in the case when ( depends on two parameters r G M and L G C, 
where r > Tq > 0 , L ^ R+: 


0 (3.12) C = Cr(L) := 7r-(L) + ldr{^ = 

We select the branch of the square root in (ITT^ to have 
ReCr(L) > 0 , ImL ■ ImCr(L) < 0 . 


2r 


For L hxed all the points Cr(L) lie in some half-plane ReCr(L) > 7 o(L) > 0, 
hence (E3D is satisfied. It is clear that 7 r(L) ~ as r —> 00 . Therefore, for 
any L ^ M there exists a constant ci = ci(L) > 1 such that 


(3.13) 


.-Ri/2| 


< \C+\^ + \C_\^ < Cir 


1 / 2 1 


Vn — -|- C-ip^ £ fbc, C — Cr-(L)- 


pr 


sa. 4 


sa.4y 


sa.4a 


sa.4b 


4. Self-adjointness: proof of Theorem 12.31 

According to the general theory of self-adjoint operators, we must show that 
the equation = LR has only the trivial solution for some (and then all) 

L G C±. To simplify our notation, we shall denote 

(4.1) C™,„(L) = C.„.„(L) = (r„,„ - L)'/2, 




If R ~ {vm,n}, then each function Vm,n can be written as 

(4.2) L) + C" ,,(p-,,(x; L)) , 

with some coefficients We have inserted the factor rUX in order that (cf. 

(EH) 

(4.3) J e ^{RG(^}. 

The matching conditions (ESI) at a: = ±1 yield an infinite system of homo¬ 
geneous linear equations for the unknown coefficients Below we set 

gL± = \/2/a±. Taking ()3.4|1 into account, we get from the condition at x = 1: 

d/4 


_L lR/2(^+ 2/r+Cm,n(L)rm,n / 

' m+l,nv'^^ -*-/ '^m+l,ri. ' . /i/- /T ^ I 

I _ ) \ 




(4.4) 


= 0- 
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sa. 5 


sa. 6 


sa. 7 


sa. 8 


sa.ll 


sa.lip 


It is convenient to multiply each equation by the factor r 


1/4 

171,71’ 


Let us also denote 


(4.5) 


+ _ ™l/2^1/4^1/4 


Qm.n =m' r. 


7n,n 771—1,7X1 


- _ „l/2„l/4^1/4 


= n' r. 


771,71 771,71 —ll 


Pm,n(L) — 


The equation (gaD and the similar equation coming from the condition (ESD 
at x = — 1 yield 


(4.6) 

+ I 2/i+Pm,n(L) 

^m+l,n^m+l,n + ,t , 

- I 2/i_pm,^(L) 

?m,n+l^m,n+l “r “ 77 (1 \ 

J_ — g 


f(~<+ p~2(/m,n(L)5 ^ + 

\^m,n ^m,n^ j ~ 'dm,n^m- 

f _ /^+ p~2(/m,n(L)'\ _j_ — ^ — 

\^m,n ^m,n^ j ' 


0 , 

0 . 


Denote by = fR(L) the inhnite matrix which corresponds to this system. In 
view of (D, we consider Tl as an operator in the space 

Removing in (gSl) the exponentially small terms, we come to a simpler sys¬ 
tem 


( 4 - 7 ) ?m+l,nCm+l,n + 2 /X+Pm,n(L)C'^ “ Oj 

(4-8) qm,n+l^m+l,n 4" + qrn,n^m,n-l ~ 0- 

Let 51' = IR'(L) stand for the matrix which corresponds to the system ()4.7j) - 
(EiHD, and also for the operator in 0 generated by this matrix. The operator 51' 
decomposes into an inhnite family of simpler operators. First of all, the equa¬ 
tions (I4.7|l (for and (14.81) (for C“„) are mutually independent. Further, 

hx any n G Mq. The equations in (gZZl) which correspond to the chosen value 
of n form a linear system in f'^(Mo) with the Jacobi matrix 

3ni^) ~ 3 ({2/i+Pm,n(L)}, {qm,n}) ■ 

In the same way, the equations in (USD, which correspond to the chosen 
value of m, form a linear system in (No) with the Jacobi matrix 

5m(L) = S ({2/i_p^,„(L)}, . 

The above reasoning shows that 

(4.9) 3i'(L) = ^®3:(L)®E*37(D. 

n m 

The original operator 51 can be written as 

(4.10) Ji(L) = Ji'(L) + N(L) 
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sa.lln 


where K = 3sf(L) is a block-diagonal matrix with 2 x 2-blocks: 
(4.11) 1V„,„(L) = '*+ 


3 




3—2(^m,7T(L) 


The last two equations elucidate the structure of the matrix fk(L). 


In the rest of the section we take L = ir G ZM. We will show that each 
term on the right-hand side of dOj) is an invertible operator in and that the 
norms of (fT)~^|| are uniformly bounded. For this purpose, we note that 

(2/i+)"^Ima+(ir) = (2p_)"^ Im (zr) = diag{pm,„(zr)}. 


sa. 13 


We have = X + iY where 


2Y^ = 


m^n 




_ 

m,n 






> 2c^ 


T 


The last inequality, with some constant c > 0, is valid for |r| > tq and for any 
m, n > 0; we have taken into account that > 'f' 0,0 = (p+ + ^-)/2 - 

By the well known estimate for the operators with sign-defined imaginary 
part, see e.g. Theorem IV.4.1 in j2], this implies that 

l|3f(»T)-‘|| < (cUm)^'. Vie No, 

and therefore 

(4.12) p'(r)-^|| =sup||afc(r)“i < . 


The norms of the blocks Nm,n{i'T) in (imi) are controlled by \pm,n{i'T)\e Ym.niir) 
and hence, are bounded uniformly in m, n G Nq. Therefore, 

P(±zr) - 3^'(±Zr)|| < C = C'(ro). 

Choosing |r| large enough, we conclude from (iroi) that 
P(±zr) -D^'(±zr)|| < P'(±zr)-'||-\ 

It follows that the operator fk(±zr) has bounded inverse, and, in particular, 
the system (HH) has only the trivial solution in 0. 

The proof of Theorem 12.31 is complete. 


5. Representation of the resolvent (Aq, - L)“^ 

5.1. Resolvent. In order to prove Theorem 12.81 we need a convenient rep¬ 
resentation for both resolvents involved in (ITTTD . Here we do this for the 
operator Aq,. We derive an analogue of the formula (6.6) (the basic formula) 
in j3] or, equivalently, (6.4) in (rUj. However, there is an important difference 
between the techniques we employ here and those in and ng, The main 
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auO 


res.lamn 


res. 3 


res. 3q 


xz 


goal in both cited papers was the direct study of the a.c. spectrum of the oper¬ 
ator corresponding to the case of one oscillator. To achieve this objective, the 
behaviour of the resolvent as the spectral parameter approaches the real axis 
was carefully studied. What we do here is to apply scattering theory, and use 
the already known results of j3] and ng. This makes our analysis much easier. 
We are able to work with the resolvents for a hxed value of the parameter L. 
We exclude L from the notation, unless to do so would be confusing. 

Let a function F E 9) have the decomposition F For any a > 0 

let us denote 

(5.1) Uqi L^ F ~ y ^ 

In the notation for v we do not reflect dependence on the parameter a. 

The operator (Aq — L)“^ can be written in an explicit form. The functions 
uo-m,n are given by the formula (EH), with ( = See ()4.1j) . It follows 

that 

(5.2) M0;m,n(±l) = (2U,n)-' f 6“^™-(t)dt; K^J(±1) = 0. 

Jr 


Each function Vm,n{x) belongs to the space Hence, the equalities 

dH are satished for it, again with = (Cm,n(F). Using these equalities and 
taking into account that [m'„.^„](±1 ) = [n(„„](±l), we hnd from the matching 
condition in (EH at the point x = 1 that 


2Cm,r 


(5.3) 


_ g — 4^771,n 

o^_i_ 


As in section Hi we let /i± 
()5.dj) yields 


^\/?7Z -|- “I" ^T?7'‘Uq.- 772—l,n ( 1)^ ■ 

Since Vrn,n the eQlintion 


(6.4) 


_ g — ‘^<>^771,71 ' ^ 

2h+Cm 


mu. 


i;m—l,n(l) 


(M0;m,n(l) “ 6 ^^"‘’"'Ii0;m,n(-1)) • 


_ g 4^777,^71 

The next step is the same normalization as in section EJ Denote 


(5-5) '^m,n ^0;m,n(il), ^m,n ^m,n ^a;m,n(il)' 

Each function Vm,n can be represented as in (14.211 . with 

We use a shortened notation for the corresponding elements in 0: 


2^ = ^ C* = {C'+U" J, 


m, n G Nq. 
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res.3r 


cs 


3e 


0 


Multiplying each equation in (j5.4jl by rm,n and writing out the similar equa¬ 
tions coming from the matching conditions at x = —1, we reduce the system 
to the form 


(5.6) 


„+ 7+ , _ -2U,ri7- 'i I „+ 7+ 

_ g— 4Cm n ^ ^m^n) l,n 

= ‘^^^+Pm,n /^+ _ -2C™.n 1 . 

1 _ e-4U.n \ rn,n ^ ) 

Qm,n+l^m,n+l + {^m,n ~ ^ 

‘^P-Pm,n _ -2Cm.n;7+ t 


_ g — 4^rn,r 


This is the non-homogeneous counterpart of the system (gSl). In order to write 
it more conveniently, we need more notation. All the operators introduced 
below depend on L and we always assnme that L ^ M. 


Dehne the operator S = S(L) : 0 by 


S : F 


rV4 
' m,n 


3 Cm,n\^ 1| 




^1/4 

' m.n 




According to (Q and (EUl), this can be written as 

(5.7) S . F \ > \^p^ n^m,n'tPrn,n^rn,n}- 

It follows from the Cauchy-Schwartz ineqnality and (ixni) that the operator S 
is bonnded. 

The diagonal operator 


T J’(h) . {Xrn,ni ^m,n} {Prn,n^rn,n^ Prn.,n^rn,n} 

acts in 0 and is nnbonnded. Its inverse is a bounded operator. 

Farther, let M = M(L) be the operator generated by the block-diagonal 
matrix, M = diag{Mm,n}, where 

(5.8) M,„,„ = (1 - e-X'"’”)-' 

Evidently, M is bonnded in 0. 

Finally, we let 

(5.9) T = T(L) : + C'm,n7’m,n)}- 

This is a bonnded operator acting from 0 into Sj. 

The system ()5.6j) can be written in the operator form 



^iZ = 2MSF, 
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3f 


matr 


reso 


auOO 


ula 


whence Z = 201 ^MSF. Here 01 = Di(L) is the operator in 0 which corresponds 
to the left-hand side of the system (15 .hp . or, equivalently, of the homogeneous 
system (14.511 . We also have X = so that 

c = Z -X = - Ip-^) §F. 

If C is found from this equation, then evidently Ua — Uq = 7C. 

Now it follows from the construction that 

(5.10) (A„ - L)"^ - (Ao - L)-^ = 7 S. 

This is the desired representation of the resolvent of the operator A^. 

5.2. On the matrix 1 P(L). It was shown in section |3] that for r large enough 
the matrix CP(ir) has a bounded inverse. This allowed us to conclude that the 
operator A^ —L has a bounded inverse for all L ^ M, and hence kerCP(L) = {0} 
for all such L. So, the operator is well-defined. However, this does not 

imply automatically that this operator is bounded in 0. We now show that this 
property is a direct consequence of the representation dsnni). Indeed, (ICTl 
implies that 

20l-^M - 7-^ = {7*7)-^7* ((A„ - L)-^ - (Aq - L)"^) 

It is easy to show that for L ^ M_|_ the operators M and T*T and §§* (acting 
in 0) have bounded inverses. This yields the desired result. 


6. Representation of the resolvent (A° - L)“^ 

Our aim here is to derive an analogue of the representation (jnni) for the 
operator A°. One possible way to proceed is to use the decompositions (Q, 
(El- However, we prefer another way, one which parallels our argument in 
sectional The calculations are easier for A° than for A^^. 

For the objects related to the operator A° we use the notation 
(6.1) u: = (A» - L)-V ~ V ~ {<,„} = 

An analogue of dEl is given by 

(6-2) <™,n(±l) = (2Cm.n)"' [ fm,n{t)dt. 

J M-|- 


Next, we derive an analogue of (Q. Taking (jSUl) and [mo.^„'](±1 ) = 0 into 
account, we get from the matching condition at x = 1: 


2Cm,r 


1 — e 
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3d 


3m 


Since ~ 'w® find, taking, as before, /i± = y/2/a±. 


2 /l+Cm, 


(6.3) 


Vm + 1< +1 (1) + --+ Vm<;rrx-l,n(l) 

’ ' ’ I _ p ^^m,n ’ ’ ’ ’ 

_ ‘^f^+Cm,n o /i\ 


_ P 2^771^71 


This is much simpler than the system D, which, of course, merely reflects 
the special structure of the operator A° as given by (EH). 

The normalization, as in section 0 reduces (EH and the similar equations 
for a: = — 1 to the form 


(6.4) 


Here 




+ , 2/i+p m,n ryo + + + ‘^f^+Pm,n 




+ 


_ P 2(^777^ 

ryo — 


L 

ryo,— ‘^f^—Pm,n T^o,— 


^0,+ 
_ P ’ 


^m.n + 




m,n+l‘^m,n+l ' i ^m,n ' ^m,n+l1 i 

’ ' ’ ' I _ p ^^m,n ’ > I » I _ p ^^m,n ’ 


WO,± ^ .^-1/4 O ; 70 ,± ^ -1/4 o /Ii'. 

m,n m,n ^0;m,n\ /’ m,n m,n ^a;m,n\ /' 

The coefficients Pm,n are the same as in (EH, being defined in (jlH). By 

EH, we have 


= {e fm,n{t)dt. 


Now we define analogues of the operators involved in the equality ()5.6|) . First 
of all, = fk°(L) is the operator in 0, defined by the infinite matrix which 
corresponds to the left-hand side of EH- The operator can be written in 
the form similar to (j4.1 Oj) : 

^ (6.5) a^°(L) = ai'(L)+X°(L) 


where 


roO 


( 6 . 6 ) 


N°(L) = diag{X^,J, 


fp+ 0 \ 

1 — (; — 2(m,n y 0 P—j 


The self-adjointness of the operator A° in Sj implies that the operator fk°(L) 
is invertible for any L ^ M. 

The operator S° = §°(L) : 0 is a bounded operator defined by 

S° : F ~ {fm,n} I—>■ {Pm,nX^'^, Pm,nX^„}; 


cf. (EH). 
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Sec 


3fc 


pr(p6ffi 


goalO 


prO.qq 


p2 


The operator M° = M°(L) is the bounded operator on 0 of block-diagonal 
form 


(6.7) M° = diag{M™,„}, Mm,n = (1 - e 


- 1 / 7 ^+ 0 
0 


Finally, let 

r = r(L) : {C+,„, C-J ^ } 

This is a bounded operator acting from 0 into S). 

As in section El we can re-write the system dEill) as 
( 6 . 8 ) (A° - L)-^ - (A° - L)-i = T° (2(3^°)-iM° - S°. 

Note that for any L 7 ^ L the operator fk°(L)“^ is bounded. The proof is the 
same as for the operator fk“^(L), see section 

7. Proof of Theorem 12.81 
7.1. The case a = 0. Here we show that 
(7.1) (A°-L)-3-(Ao-L)-3e6i. 

Recall that in the notation of (EU and (10 

(Ao - L)-^F = {no,™,4, (A° - L)-'^ = 

where the functions no,m,n,Momn given by the equations dO and (ITTUD 
respectively, with ( = Cm,n- So, both operators are diagonal. Denote by 
fheir components, and let 

Q = diag{g™,4 =(A° - L)-1 - (Ao - L)-\ 

‘ O - $ 

^m,n ^m,n- 

According to (ITTTl) . each Qm,n is a rank one operator: 


(7.3) Qm,n ■ ^ ^ (2Cm,n) 9m, ni^') I 9m,ni^^ f m,ni^^dt 

jR 

where gm,n{,x) = It follows from ()7.3jl that 

||g™,n|| = (2|C™,n|Re W)-' < C'r-;„, C = C'(L). 

The norms of and ^ ^e easily estimated (actually, ||d’m,n|| can 
be calculated explicitly, since this is a convolution operator). By the “Schur 
test”, the norm ||K|| of an integral operator in with the kernel K{x,t) can 
be estimated as 
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s 


end. 0 


end. 1 


end 


Applying this to the operators (uni) and innui) . we hnd that 


Id) II lld)° II < 


Furthermore, the components of the operator (Aq — L) ^ — (Aq — L) ^ are 

^m,nQm,n + ^m,nQm,n^rn,n Q'm,ni^rn,n) ' 

The norm of this operator does not exceed and, since its rank is no 

greater than 3, its trace class norm does not exceed By ()2.5j) . these 

numbers form a convergent double series, and hence, (EU is established. 

Note that the exponent 3 in (EH) can not be replaced by 2. 


7.2. Difference between the right-hand sides in (Ib.KIj) . (Ib.Hj) . To shorten 
our notation, let us denote 

(7.4) H = T (2ai-iM - S, H° = T° ((23^°)-^M° - S°. 

Here we show that 

(7.5) T;=H°-H=g6i. 

Since all the inverse operators appearing in (El are bounded, we only need 
to check that 


7 ° - T, S° - §, - M, 0^° - G ©1. 

Here each operator has a block-diagonal structure, with (2 x 2)-blocks, and it 
is sufficient to estimate the operator norm of each block and to verify that the 
corresponding series converge. 

For the operators T° — T and §° — S the result immediately follows from the 
dehnitions of the operators involved and the estimate (El- For the operator 
M° — M the result is evident from the comparison of (El and (El- Finally, 
for — fk the result is implied by (14.1011 and (jti.5|) . if we take into account 
evident estimates of the norms of blocks Nm,n in (I4.1H1 and in ( 10 . 0|) . 

7.3. End of the proof. Unfortunately, the desired inclusion (imi) is not 
implied by mil and mi automatically and we need an extra argument in 
order to dualize the proof. 

Let us denote G = (Aq — L)“^. Using also the notation Q, T as in (El, 
El, we can re-write the equalities (ICTIll and (El as 

(A„-L)-^ = G + H, (A° - L)-i = G -I- H -I- Q -f T. 
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We already know that \1/ G 6i and (G + Q)^ — G 6i. Therefore, the 
following equality is satished modulo a trace class correction: 


(A° - L)-' - (A„ - L)-i = (G + Q + H + T)3 - (G + H)3 
= (G + Q + H)3-(G + H)3 ( mod6i) 

= ((G + Q)2 - G^) H + (G + Q)H(G + Q) - GHG 
+H ((G + Q)2 - G") + QH^ + HQH + H^Q ( mod 6i). 


Removing the parentheses, we come to the sum where each term involves one 
of the products QH, HQ, QGH, HGQ. Taking into account the structure of 
the operator H, we see that it is sufficient for us to prove that the operators 


QT, SQ, Q(Ao-L)-iT, §(Ao - L)-iQ 


(7.6) 


are trace class. All these operators have block-diagonal form, with the blocks 
given by explicit formulae implied by the corresponding dehnitions. For in¬ 
stance, according to (EH) the operator QT transforms the number sequence 
into the sequence of functions {wm,n} where 




An elementary calculation shows that the integral here is of order 0(6“'*'"*’"). 
This happens because the function gm,n(t), see (|7.d|l . is concentrated around 
the point f = 0, while is concentrated around t = ±1, and all the 

three functions decay exponentially when t moves away from the corresponding 
center. Clearly, this estimate implies that QT G ©i. The proofs for the other 
operators in dZH are similar. 

The proof of is complete. 
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